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Detailed analysis of spectral line broadening and variations in relative intensities of hyperfine 
spectral components due to optical pumping is presented. Hyperfine levels of sodium 3pi/2 and 3^3/2 
levels are selectively excited in a supersonic beam at various laser intensities under the conditions 
when optical pumping time is shorter than transit time of atoms through the laser beam. The 
excitation spectra exhibit significant line broadening at laser intensities well below the saturation 
intensity, and redistribution of intensities of hyperfine spectral components is observed, which in 
some cases is contradicting with intuitive expectations. Theoretical analysis of the dynamics of 
optical pumping shows that spectral line broadening depends sensitively on branching coefficient 
of the laser-driven transition. Analytical expressions for branching ratio dependent critical Rabi 
frequency and critical laser intensity are derived, which give the threshold for onset of noticeable 
line broadening by optical pumping. The critical laser intensity has its smallest value for transitions 
with branching coefficient equal to 0.5, and it can be much smaller than the saturation intensity. 
Transitions with larger and smaller branching coefficients are relatively less affected. The theoretical 
excitation spectra were calculated numerically by solving density matrix equations of motion using 
the split propagation technique, and they well reproduce the observed effects of line broadening and 
peak intensity variations. The calculations also show that presence of dark (i.e., not laser- coupled) 
Zeeeman sublevels in the lower state results in effective branching coefficients which vary with laser 
intensity and differ from those implied by the sum rules, and this can lead to peculiar changes in 
peak ratios of hyperfine components of the spectra. 

PACS numbers: 32.70.-n, 32.70.Jz, 32.80.Xx 



I. INTRODUCTION 

Optical pumping is a well known phenomenon, which 
is usually associated with redistribution of population 
within hyperfine (HF) components or Zeeman sublevels 
of the ground state due to coupling by resonant light 
fields Optical pumping is being exploited in var- 
ious applications, like cooling below the Doppler limit 
@> HI j vibrational excitation of molecules in the electronic 
ground state [3| , orientation and alignment of atomic and 
molecular ground states [f|, etc. When optical pumping 
is involved in the control of quantum states, it is usu- 
ally associated with large laser intensities exceeding the 
saturation limit 0. Therefore, the populations of quan- 
tum states depend nonlinearly on laser intensities and 
the excitation spectra are affected by power broadening 

In present study we are concerned with lineshape ef- 
fects due to optical pumping in the weak excitation limit. 
Specifically, we measure laser excitation spectra of the 
3pi/2 and states of Na in a supersonic beam. Cou- 

pling of the F" — 1 and F" — 2 levels of the ground state 
with different HF components of the upper states allows 
us to study two-level systems with different branching 



Nozzle Skimmer 1 

d = 0.4 mm A = 2.0 mm 



Skimmer 2 
d> = 1 .5 mm 



Aperture 
b = 3 mm 



*Fock Institute of Physics, St. Petersburg State University, 198904 
St. Petersburg, RUSSIA 

tlnstitute of Semiconductor Physics SB RAS, 630090, Novosibirsk, 
RUSSIA 



1_J 



21 23 25 
— I 1 1 



L, cm 

FIG. 1: Collimation of the sodium beam by skimmers and 
apertures. Laser beam crosses the atomic beam at right angles 
25 cm downstream from the nozzle. 



coefficients. The smaller the branching coefficient, the 
more population irreversibly leaves the two-level system, 
and vice versa. At very low laser intensities the exci- 
tation spectra do not reveal any abnormalities. When 
laser intensity is increased but still below the saturation 
intensity, essential modification of the excitation spectra 
is observed: most of the hyperfine spectral components 
exhibit additional broadening while their intensity ratios 
cease to obey the line strengths rules. Note, that usually 
line broadening is considered to be a strong-field effect 
due to power broadening at laser intensities above the 
saturation intensity 0, El] • 

Since the experiments were performed at low number 
densities of sodium atoms (n^s ~10 10 cm -3 ), lineshape 
modifications by radiation trapping can be disregarded 
0, [1] • We attribute the observed lineshape effects to op- 
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tical pumping, which leads to depletion broadening of 
spectral lines [3, [To| . If transit time r tr of atoms through 
the laser beam is much larger than lifetime T na t of the 
3p state, populations of levels and the associated fluo- 
rescence signals can become nonlinear on laser intensity 
lias long before the saturation limit is reached (i.e., at 
lias <C I sat)- Due to interaction with laser field the 
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where T nat = 2tt ■ Av nat — 1/ T nat is the natural width (in 
units of angular frequency [sec -1 ]) of the excited state e, 
S is the laser detuning from the line center, and Q is the 
Rabi frequency of the transition. The width T" is equal 
to the rate of photons spontaneously emitted from state 
e. If the level system is partially open, only the fraction II 
of the spontaneous transitions will return the population 
to the initial state g; the fraction 1 — IT associated with 
decay to levels other than g will be lost from the (g, e)- 
system during each excitation-emission cycle. The rate of 
such pumping is obviously T pump = (1 — II) T" . Hence, 
the pumping time can be written as 
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If the transit time r tr is long, such that r tr > Tpump — 
Tpump(S — 0), the population of the (g, e)-system will 
be fully depleted during interaction with the laser field. 
In terms of Rabi frequencies the condition for pop- 
ulation depletion can be rewritten as f2 > fl cr = 
QsatyJ^Tnat/ {r t r (1 — II)), where saturation Rabi fre- 
quency is Q S at = Tnat/v^ [H . Note, that the parameter 
1 / \jTnatTtn which was considered in [f| [l2j as the param- 
eter associated with saturation due to optical pumping in 
the case of open level systems (i.e., no population return 
form state e to state g), is identical to our critical Rabi 
frequency Sl cr in the limiting case of 11=0. 

If the weak excitation limit is combined with long in- 
teraction times of atoms with the laser field, such that 
T tr >> T nat , the value of critical Rabi frequency is small 
(fl cr < Qsat) and broadening and saturation of spectral 
lines can be observed at laser intensities well below the 
saturation limit, long before power broadening starts af- 
fecting the lineshapes. 



II. EXPERIMENT 

The experiment was performed in a supersonic beam 
of Na atoms (see Fig. 1). Two skimmers and an entrance 
aperture of the excitation zone collimate the beam with 
flow velocity !/j to a small divergence angle thus re- 
ducing the Doppler width for excitation perpendicular to 
the beam axis to Av D = t>/$/(2A), where A is the laser 



FIG. 2: Hyperfine energy levels of the 3s and 3p states of Na. 



wavelength. The divergence angle $ was varied between 
0.67° and 0.92° by using the entrance aperture b of either 
2 mm or 3 mm dia. The laser beam crosses the atomic 
beam at right angles, and it is linearly polarized parallel 
to the molecular beam axis z, which is also the quantiza- 
tion axis. Only Zeeman sublevels with identical quantum 
numbers mp are coupled by the laser field due to the se- 
lection rule ArriF = 0, while for transitions between lev- 
els with the same F the transition mpn = <-> mpr = 
is forbidden. The number density of atoms in the beam 
was chosen sufficiently low (< 10 10 cm -3 ), thus ensuring 
that the beam is not optically thick and effects of radia- 
tion trapping and photon reabsorption Q can be safely 
neglected. An important consequence of optical trans- 
parency of the beam is that absorption P(Av L ) (the 
total number of photons absorbed per second) and ex- 
citation J (Av L ) (the integrated over frequencies flux of 
emitted photons in the direction of observation) profiles 
as function of the laser detuning Avl do not vary in the 
interaction volume defined by the crossing atomic and 
laser beams. Both profiles are proportional to the in- 
tegral (over the interaction volume and HF sublevels) 
population of the excited state. 

The 3s — > 3p transition was excited using a single mode 
cw radiation source (Coherent CR-699-21 dye laser) with 
linewidth of 1MHz. The fluorescence emitted by Na 
atoms was collected into two fiber bundles at the angles 
of 90° and 45° with respect to the directions of the axis of 
the molecular beam, laser beam, and laser polarization. 
The fluorescence light was guided via the fiber bundles 
to two photomultipliers, and the signals proportional to 
J (Av L ) were registered using photon counters. The re- 
sulting excitation spectra of the 3p state were recorded as 
a function of laser detuning Av^ . The arrangement with 
two different simultaneous detection geometries allowed 
us to verify that radiation trapping, which is strongly 
anisotropic with respect to the direction of the observa- 
tion, docs not affect the measured spectra. It also al- 
lowed us to rule out the influence of polarization effects 
on variations in lineshapes and relative line intensities. 

The mean flow velocity Vf of atoms in the beam was 
measured to be 1160 m/s. For excitation perpendicu- 
lar to the atomic beam axis the apertures of the excita- 
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FIG. 3: Excitation spectra of the 3s 1/2 , F" = 2 ^ 3p 1/2 ,F' = 
1,2 transitions in Na. Residual Doppler width due to finite 
collimation angle is Az/£>=11.2 MHz (at b — 2 mm). The 
expected peak ratio is 1:1. Saturation intensities of the lhs 
and rhs components are 7.5 and 12.5 mW/cm 2 , respectively. 



tion zone of b = 2 and 3 mm correspond to the residual 
Doppler width Av D = 11.2 and 15.9 MHz (FWHM), re- 
spectively. These should be compared to the excitation 
perpendicular to the natural width of Ai/ nat of 9.8 MHz 
(r no t=16.23 ns for the 3p 3/2 [3). 

Figure 2 shows the hypcrfinc energy levels of the 3s and 
3p states. The excitation spectra were obtained by scan- 
ning the laser frequency across the 3si/ 2 — > 3pi/2 an d 
3^1/2 — ¥ 3^3/2 transitions. The HF splittings are larger 
than both the Doppler width and the natural width for 
all but one pair of components {ipz/2 F' = and F' = 1). 
The measurements for the Di-line (A=589.593 nm) were 
performed with the aperture b—2 mm. In the case of 
the D2-line (A=588.996 nm), 6=3 mm was used. Ra- 
dius of the laser beam was r; as =1.5 mm, which corre- 
sponds to the transit time Tt r — 2 x ri as /vf=2.65 ps at 
t>/ = 1160 m/s. Thus, the transit time is by more than 
two orders of magnitude larger than the natural lifetime 
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FIG. 4: Excitation spectra of the 3s 1/2 , F" = 2 -> 3p 3/2 ,F' = 
1, 2, 3 transitions in Na. Residual Doppler width due to finite 
collimation angle is A^_d = 15.9 MHz (at b = 3 mm). The ex- 
pected peak ratio is 1:5:14. Saturation intensities of the three 
components are 37.4, 12.5, and 6.2 mW/cm 2 , respectively. 
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FIG. 5: Excitation spectra of the 3si/2, F" = 1 — * ?>Pi/2, F' = 
1,2 transitions in Na. Residual Doppler width due to finite 
collimation angle is Ai/d = 11.2 MHz (at b — 2 mm). The 
expected peak ratio is 1:5. Saturation intensities of the lhs 
and rhs components are 37.4 and 12.5 mW/cm 2 , respectively. 



of the 3p state. 



III. SPECTRA 

Figure 3 shows the measured excitation spectra of the 
3pi/2 state from the F"—2 sublevel of the ground state 
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FIG. 6: Line strengths (square frames) and branching 
ratios ILj (circular frames) for (a) 3si/2 — * 3pi/2 and (b) 
3si/2 —* 3p 3 / 2 hyperfine transitions. 



at various laser intensities. The spectra exhibit two 
peaks corresponding to the excitation of the 3s]/2, F" = 
2 — > 3pi/2,F' — 1,2 HF transitions. The spectra of 
Fig. 3 were measured at the divergence angle of the 
atomic beam of i? = 0.67°, which corresponds to resid- 
ual Doppler width of Av D = 11.2 MHz (see Sec. VI for 
details on Dopller lincshape). The spectrum of Fig. 3(a) 
was measured at a very low laser intensity of 21/iW/cm 2 . 
Both HF components appear equally strong, which is 
obviously due to equal line strengths of both HF tran- 
sitions The lineshapes are determined by a com- 
bined effect of natural (Av nat — 9.8 MHz) and Doppler 
broadening. In Fig. 3(b) the laser intensity has been in- 
creased by a factor of about 70 compared to Fig. 3(a) 
to the value of 1 mW/cm 2 , which is still much smaller 
than the saturation intensity of both HF components 
(7.5 and 12.5 mW/cm 2 ). One can observe that the 
F" = 2 — * F' = 2 (rhs) component has become some- 
what smaller than the F" = 2 — ► F' = 1 (lhs) com- 
ponent. Peculiarly, when the laser intensity is further 
increased to 25 mW/cm 2 [Fig. 3(c)], the rhs component 
becomes somewhat larger than the lhs component, while 
the the widths of the peaks (Au = 75 MHz) are substan- 
tially larger than the width Av sat = 16.5 MHz expected 
from saturation broadening at this laser intensity. 

The 3s 1/2 ,F" = 2 -► 3p 3/2 ,F' = 1,2,3 excitation 
spectra are shown in Fig. 4. The relative peak inten- 
sities match the theoretical line strengths of individ- 
ual HF transitions when laser intensity is very small 
[15/iW/cm 2 , Fig. 4(a)]. When laser intensity is increased 
to 1.5/iW/cm 2 , which is still below the saturation in- 
tensity, the relative intensities of the components corre- 
sponding to the excitation of the F' = 1 and F' = 2 
HF levels are smaller than expected from the theoretical 
line strengths [Fig. 4(b)]. When laser intensity is close to 
saturation intensity, the F' = 1 and F' — 2 peaks are so 
weak compared to the F' = 3 peak that it is ambiguous 
to attempt analysis of their linewidth. 

The 3si/2,F" = 1 — * 3pi/2,F' = 1,2 excitation spec- 
tra are shown in Fig. 5. Like in the case of Fig. 3, also 
here a significant broadening is observed at laser intensi- 
ties below the saturation limit. In contrast to Fig. 3(b), 



however, the lhs peak corresponding to the excitation of 
the F'=l component of the upper state grows monoton- 
ically as compared to the rhs peak. 



IV. THEORETICAL LINE STRENGTHS AND 
SATURATION INTENSITY 

At very small laser intensities the strengths of individ- 
ual peaks in the excitation spectra shown in Figs. 3-5 cor- 
respond to the respective theoretical line strengths 
of individual HF transitions i = {F" — * F'} within the 
Di (j = 1/2) or D 2 (j = 3/2) lines [14]. The values of 

Si are directly related to the reduced matrix elements 
of the transitions and the partial natural width of the 
respective transition: 



r (i) = 4uf 1 aU) 
nat 3Hc 3 2F' + l 1 



Sf = \{l/2,F"\\D\\j,F') 

(3) 

Intensity of the component i is proportional to its line 
strength because the product T^ at (2F' + 1) regulates the 
photon flux of this component under the conditions of 
thermodynamic equilibrium ■ Figure 6 shows the the- 
oretical line strengths (square frames) in units of 
the reduced matrix element | (3s||D||3p) | 2 of unresolved 



3s -> 3p transition, i.e., = >\ (3s\\D\\3p) | 2 . The 
values of S^' are normalized such that Hj = 8 [l4| . 
Thus, the theoretical values of peak ratios can be directly 
taken from Fig. 6, and they agree with the experimental 
observations at very small laser intensities [see Figs. 3(a), 
4(a), and 5]. 

Saturation intensity of each hyperfine transition de- 
pends on the natural width of the transition Av nat and 
the branching ratio Q (see also in Sec. V): 
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The values of the hyperfine branching coefficients (cir- 
cular frames in Fig. 6) are easily obtained from the re- 

duced line strengths S\ . Note, that all HF transitions 
have the same natural width of Ai/ raat =9.8 MHz. 

The saturation intensity given by Eq. gives the 
limiting laser intensity after which stimulated transitions 
start transforming the excitation spectra Q. The low- 
est laser intensities used in measurements of the spectra 
shown in Figs. 3-5 are much smaller than saturation in- 
tensity of any of the HF transitions. Therefore one 
naturally expects the peak ratios to be in accordance 
with the line strengths of Fig. 6 and line width to cor- 
respond to residual Doppler width determined by beam 
divergence. This agrees with the observations made for 
the smallest laser intensities. 

When laser intensity is increased by a factor of about 
100, it is still well below the saturation intensity I, 
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Nevertheless, a curious transformation of the spectra is 
observed: widths of the peaks increase, and peak ratios of 
the HF components change. Interestingly, the strongest 
F" = 2 -> F' = 3 component of the D 2 -line (Fig. 4) is 
not affected by broadening at all although its saturation 
intensity 7^=6.2 mW/cm 2 is the smallest. Intuitively, 
one would expect the transition F" — 2 — > F' = 3 to 
be the first one that is affected by broadening when laser 
intensity is increased. 

Not only the widths are affected. Relative intensities 
of the HF peaks change as laser intensity is increased. 
At the first glance it seems that relative intensities of 
components with smaller branching ratios H should de- 
crease when optical pumping becomes non-negligible, as 
less population returns to the lower laser coupled level 
than it does for levels with larger 11^. This is clearly the 
case in Figs. 3(b) and 4, but not in the case of Fig. 5. 
Moreover, Fig. 3 shows another unexpected feature: af- 
ter the intensity of the peak with smaller H has initially 
decreased with respect to the peak with larger 11^ [cf. 
Figs. 3(a) and 3(b)], a further increase of laser inten- 
sity leads to increase of the peak with smaller H [cf. 
Figs. 3(b) and 3(c)]. Explanation of these observations 
requires a detailed analysis of the dynamics of optical 
pumping. 



V. DYNAMICS OF OPTICAL PUMPING AND 
ITS EFFECT ON THE FLUORESCENCE 
SIGNALS 

The measured fluorescence signals are affected by var- 
ious factors, like the detection efficiency and geometry. 
The spectral components are excited and detected at very 
close wavelengths under identical conditions, therefore it 
can safely assumed that the detection efficiency is equal 
for all of them. Since we are interested in relative in- 
tensities and widths of the components, it is sufficient 
to consider the fluorescence signals that are proportional 
to the total number of photons emitted by atoms at all 
times in all directions. 

We consider the following model problem. Two-level 
atoms with the ground state g and the excited state e 
propagate along the z-axis with the flow velocity of the 
beam vt. The atoms cross the laser beam with radius 
ri asi frequency to and Gaussian intensity distribution 



I (z) = 7, as exp (-z 2 /rf as ) ; r tr = 2r las /v f . (5) 

It corresponds to Gaussian switching of Rabi frequency 
il of the g — e transition: 



fl (t) = n exp (-2t 2 /"4) ; t = z/v { . 



(6) 



The value fi = E (g\d z \e) is the Rabi frequency of the 
g — e coupling in the center of the laser beam, which is 
linearly polarized parallel to z-axis. 



A. Evaluation of the fluorescence signal 

In what follows we shall assume that transit time is 
much larger than lifetime of the upper state, r tr >> r nat , 
which is true for the parameters of our experiment. This 
allows us to use the adiabatic elimination for the non- 
diagonal density matrix element p eg [l5| : 



Peg (t) 



in (i) 



(n g (t) - n e (i)) 



(7) 



The above equation relates p eg (t) to the populations 
n g(t) — Pgg(t) an d n e (t) — p ee (t). The decay rate 
T e = l/r.nat gives the natural width of level e, while 
5 = 2-k/S.vl is the laser detuning. With p eg defined by 
Eq. , the time evolution of the populations is given by 
simple balance equations: 



dt 
d 

dt n ° 



T e n e + r (t) (n g 
r (t) (n e 



nr e n e 



n g ). 



(8) 
(9) 



The first equation describes the population loss from level 
e via two processes: (i) spontaneous decay at the rate 
r e , and (ii) stimulated emission at the rate equal to the 
optical pumping rate r (t): 



r(i) 



n 2 jt) 
'WTf 2 , 



(10) 



The population of level g is affected by three competing 
processes: (i) photon absorption at the rate r (t)) result- 
ing in the population of level e, (ii) return of population 
from level e to level g due to stimulated emission, and (iii) 
return of population from level e to level g due to spon- 
taneous emission. The rate of the latter is determined 
by the branching coefficient II of the given HF transi- 
tion. The branching coefficients are normalized such that 
n=0 for an entirely open system (no spontaneous return 
from level e to level g) and n=l for a closed system (no 
transitions outside outside the g — e system). 

The initial conditions of Eqs. © and © follow from 
the requirement that initially all the population is in 
level g while level e is not populated: n g (t = — oo) = 1; 
n e (t = — oo) = 0. The assumption r tr > T nat leads to 
a further simplification of Eqs. © and |SJ in the weak 
excitation limit, when r (t) < 0.5r e . As weak excitation 
we understand excitation at laser intensities smaller than 
the saturation intensity given by Eq. i.e., when Rabi 
frequency of the transition does not exceed the saturated 
value, Qo < &sat, where Q S at = r e /y2 [3j. In that case, 
the adiabatic elimination implies that dn e /dt=0 [la ], and 
Eq. © immediately yields 



r e n e + r (t) (n g - n e ) = => 
r (t) 

=> n e (t) - -y- (n g (t) - n e (t)) 



(11) 
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Equation ([9]) can then be transformed into the form 



= -r(t)(l-H) n_ (t); 
n_ (t) = n g (t) - n e (i) . (12) 



The above equation can be comfortably used for the eval- 
uation of the fluorescence signal J. Integration of both 
sides of Eq. © yields the total number of spontaneous 
photons emitted by the excited atoms: 



which is given by the ratio of transit time Tt r and pump- 
ing time Tpump- The latter was already discussed in 
Sect. I [Eq. |J2J)], and it has the meaning of optical pump- 
ing time at resonant excitation (<5=0). Importantly, 
the parameter P pump can be large even at laser inten- 

3> 1 when 



sities well below the saturation limit: P v 
T tr > Tm/m and fi < tt sat . 



B. Line broadening by optical pumping 



J = r„ 



dtn e (t) = / dtr (t) 



(13) 



When the pumping parameter is small (P pU mp <C 1), 
equation (|17[) simplifies to yield the ordinary Lorentz line- 
shapes: 



Integration of Eq. 
(pf yields 



T2l) and combination of the result with 



Jl (S) = R 
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(18) 



(1 -II) J = 1 - n- (t = oo) = 1 -n g {t = oo). (14) 

The above expression has a straightforward physical 
meaning: the number of spontaneously emitted photons 
on transitions outside the g — e system is equal to the 
total loss of ground state population during interaction 
with the laser field. 

Using Eqs. (fT2| and (fl"4f . we can derive the fluores- 
cence signal in an explicit analytical form: 
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Since we consider the case of weak excitation, when 
r (t) < 0.5r e , the integral R in Eq. JT5|) further simplifies 
to the form 



R = 



dtr (t) = 



V^FeTtr 



O 2 
"0 



2 AS 2 + T 2 e ' 

n < n sat = r e /V2. 



(16) 



Dependence of the fluorescence signal on the laser detun- 
ing S = can now be rewritten as 



When Pp Um p is increased, Eq. (fT8|) no longer holds and 
Eq. (fT7|) must be used. An almost 10-fold increase of the 
linewidth is observed as P pump is increased form 0.1 to 50 
(see Fig. 7). Such broadening has a simple explanation. 
Consider a near resonant case, when Av^ = 8/2ir w 0. 
Starting from values P pU mp ~ 1 the atoms spend suf- 
ficient time in the laser filed for the population of the 
ground state to be depleted, n g (t — oo) ~ 0. Deple- 
tion of level g is associated with the emission of a fixed 
number of photons 1/(1 — II) (see Eq. ([14]) ). Hence, op- 
tical pumping saturates the observed signal / (5 « 0) via 
depletion saturation, provided that T pump (5) > r tr . Fur- 
ther increase of P pU mp cannot increase the number of 
photons emitted upon excitation at the line center. At 
the same time, the number of photons emitted upon ex- 
citation in the wings continues increasing with P pum p un- 
til depletion saturation is reached at consecutively larger 
laser detunings 6. Therefore linewidths in the excitation 
spectra will increase with P pU mp, and lineshapes will ex- 
hibit the characteristic flat-top peaks at large Pp Um p- 

The relative increase of the width Avop — Sop/2tt 
(FWHM) of the line profile affected by optical pumping 
as compared to the natural width can be easily obtained 
fromEq. {TTJ: 
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The above equation shows that excitation spectrum 
strongly depends on the pumping parameter Ppump ? 



Variation of the width with P pU mp is shown on Fig. 8. As 
can be seen, the broadening becomes noticeable at about 



p 



1, and further increase of the width scales as 



square root of P pU mp for large values of P pU mp- 

The condition P pum p > 1 for broadening by optical 
pumping can be reformulated in terms of Rabi frequen- 
cies, i.e., Rabi frequency of the laser-driven transition 
must be larger than some critical value f2 cr : 



7 



SI > fi r 



i7 cr — ^saM 



2r„ 



rtr (1 - n) 



(20) 
(21) 



Inequality (|2"0"1) generalizes the results obtained in 0, [l2j 
for the limit of entirely opened level systems with II = 0, 
which is often used as an approximation in the case 
of molecules with many possible rovibronic transitions. 
Broadening by optical pumping turns out to be sensi- 
tively dependent of the branching ratio II. In the limit 
of closed level system (IT = 1) fi cr is formally equal to 
infinity. It is not surprising: if there are two isolated 
quantum states, then no pumping can occur regardless 
how strong is the exciting laser field. 

It is important to note that the above described broad- 
ening mechanism differs from the classical textbook ex- 
amples of saturation and power broadening Q . In [|| , the 
saturation broadening is attributed to the strong field ef- 
fects, when light-induced pumping rate becomes compa- 
rable to the relaxation rates, while the power broadening 
is attributed to a considerable Rabi flopping frequency. 
The saturation parameter in both cases is given by the ra- 
tio of pumping rate to relaxation rate. In our considered 
case, in contrast, the pumping rate is very small, laser 
intensity is well below the traditional saturation inten- 
sity given by Eq. (TJ|, yet a notable depletion of level g is 
reached via optical pumping in a partially open two-level 
system due to long interaction time with the laser field; 
the line broadening thus occurs in the weak excitation 
limit, when n e (t) <C n g (t). 

Line broadening by optical pumping is thus the dom- 
inating line broadening mechanism when Q C rit < SI < 
Cl sa t- At > Cl sa t the rate of laser-induced transitions 
exceeds the spontaneous transition rate. When the split- 
ting (~ SI) of laser dressed states exceeds their widths, 
the population is equally shared between the levels e and 
g [2|, |l6[ . The pumping time T pump (0) then stabilizes 
at 2/ (r e (1 — II)) and becomes independent on further 
increase of laser intensity. 




-10 -8 -6 -4 -2 2 4 6 8 10 

Reduced detuning 

FIG. 7: Signal J (Eq. (|17[) ) as a function of the reduced de- 
tuning <5/r e for different values of the pumping parameter 
Ppump (shown as labels of the curves) . 
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FIG. 8: The reduced width Au pU m / Au„ at given by Eq. (fl9 
as a function of the pumping parameter P pU mp- 



have replaced the natural width T e by the partial nat- 
ural width r e II of the given transition. This can be 
done because Rabi frequency SI = Ed/H involves the 
dipole element d = (g\d z \e) associated with the par- 
tial natural broadening r e il = 4uj 3 d 2 /3hc 3 Since 



lias = E c/8n, we obtain Eq. 
equation (|20p can be rewritten as 



Using Eq. 



C. Critical laser intensity 

It is useful to rewrite Eq. (|20[) in terms of laser inten- 
sities, since those are usually used by experimentalists. 
We shall therefore express the laser intensity in terms of 
the transition Rabi frequency as follows: 



1-r.r — 



in 2 he 



4r„ 



V^3A3 Ttj , n (i _ n) ^T tr (l-uy 



(23) 



Simultaneously, the pumping parameter can now be 
rewritten in terms of ratio of actual laser intensity to 
critical laser intensity: 



ATT 2 hc ft 2 

3 a 3 fin' 



(22) 



The above equation is well known for closed systems with 
11=1 [3]. In the case of a partially open systems we 



Ppum — Ilas/Ic 



(24) 



Equation (|23|) shows the relation between the critical 
laser intensity and the traditional saturation intensity 
[Eq. Q], which gives the limit for onset of power broad- 
ening. 
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Note, that formally I cr — > oo when II — > (completely 
open system). The limit of II — ► corresponds to d — > 0, 
i.e., to forbidden optical transitions. Therefore, even a 
very small transfer of population from level g to level e 
will require extremely large laser intensity. 



VI. RESIDUAL DOPPLER BROADENING 

Besides the broadening due to optical pumping, the 
spectral lines are also affected by a small but non- 
negligible Doppler broadening due to finite divergence 
angle § of the atomic beam. Such divergence is associ- 
ated with non-zero velocity components in the direction 
of the laser beam for atoms moving not exactly paral- 
lel to the atomic beam axis. For atoms experiencing 
the Doppler shift Aw the laser detuning 8 will transform 
into the detuning 6 + Aw. A corresponding replacement 
5 — > 6 + Aw should therefore be done in Eq. (fl7|) . The re- 
sultant profile of a line in the excitation spectrum is thus 
given by a sum of profiles (|17p resulting from absorption 
of laser photons by atoms of different velocity groups. 
If the probability of atoms to have a velocity leading to 
the Doppler shift Aw is Pd (Aw), then the resultant line 
profile is given by the integral 



Jres 



dAwP D (Aw) J (6 + Aw) . 



(25) 



The analytical form of the function Pjj (Aw) for effusive 



beams has been derived in [17[ , while for the case of su- 
personic beams it will be analyzed in detail in 18]. This 
analysis builds on the following assumptions: (i) nozzle 
diameter d is small compared to the diameter b of the 
entrance aperture and distance L from the nozzle to the 
excitation zone (see Fig. 1); (ii) divergence angle z? of 
the atomic beam is small; (iii) size of the excitation zone 
~ ri as is small compared to the distance L, and the distri- 
bution of atoms within ~ ri as is uniform; (iv) the velocity 
distribution in the direction perpendicular to the atomic 
beam axis is due to the divergence of the beam with the 
axial velocity distribution F(v). The distribution func- 
tions F(v) for various kinds of beams can be found in 

m- 

Leaving out the somewhat lengthy detailed derivation 
of the function Pjj (Aw) to the forthcoming paper [l8| . 
we shall give here the final form of the most essential 
core (\Ai>\ < Az/p/1.5) part of the distribution function 
for the supersonic beam: 



The values of the parameters U/, A, b, d, and L are given 
in Sect. II. Note, that the function ([26]) deviates strongly 
from the Gaussian function, which is usually associated 
with Doppler profiles. The frequency dependence of Pd 
in the wings of the spectral line (\Ai/\ > Avu/l.b) differs 
from that given by Eq. (|26|) . Nevertheless, in our case it 
is sufficient to use only the core part of Pjj. Since Avd 
is comparable with Av natl the natural broadening out- 
competes the exponentially small wings of the Doppler 
profile at large Av [Hi ]. 



VII. RESULTS AND DISCUSSION 

Calculations of the theoretical spectra are performed 
in two steps: (i) solution of the evolution problem for an 
individual atom excited by linearly polarized laser field 
detuned by Av = 6/2tv, and (ii) calculation of the resul- 
tant line profile by performing the convolution (|25p . The 
first step is performed by modeling quantum dynamics 
of individual pairs of Zeeman sublevels mp within the 
F"m,F" — * F'mpi HF transition during coupling of the 
levels by the electrical field of laser light distributed as 
| E |= i? cxp (— z 2 /2rf as ) | e 2 cos(wi) |. Correspond- 
ingly, the spatial distribution of Rabi frequencies of indi- 
vidual HF transitions also follow the Gaussian distribu- 
tion: 



n tm ' = £ exp (-z 2 /2rf as ) {F"m FI ,\d z \F'm F ,) . (28) 

It is convenient to introduce the reduced Rabi frequency 
fl re d associated with unresolved 3s — 3p transition: 



E 



n red = -£y/\(3 S \\D\\3 P )\ 



(29) 



Rabi frequencies of individual Zeeman components can 
then be calculated from Q re d using the known line 
strengths S\ given in Fig. 6 and the 6j symbols 



n 



(m) 



F" 1 F' 
—itif" mpi 



(30) 



where indexes i and j stay for the chosen HF component 
and the chosen j level of the upper state, respectively. 

The values of £l r ed used in the calculations may be 
obtained from their relation to the laser intensity Ii as 
[see Eq. 



with 



p( cor) (A,) = -^Jl-A^/A^, (26) j^^Vk, (31) 

where g% p = 3 is the statistical weight of the 3p state. 
Note, that in the calculations of theoretical spectra we 
Vf d b + d used £l re d as the only fitting parameter. The theoretical 

values of J; as given in Figs. 3-5 were calculated from the 
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fitted n re d values using Eq. (|3~Tj) . and they are in a good 
agreement with the experimental values calculated from 
measured laser power and radius of the laser beam. 

For a qualitative interpretation of the experimental re- 
sults it is helpful to consider a simplified model, in which 
the populations of Zeeman components evolve indepen- 
dently, and the resulting signal J is simply a sum of in- 
dividual signals J m with the same mpn = mpi = m: 
J(Av) = ^J m (Af). In the reality, however, Zeeman 

m 

sublevels are subject to the spontaneous emission on 
transitions with Amp = ±1. As will be shown below, 
such cascading can significantly change the excitation 
spectrum as compared to the simplified treatment when 
couplings with Amp = ±1 are neglected. 

In order to account for cascading, we have elaborated 
an accurate numerical algorithm allowing the inte grat ion 
of equations of motion for the density matrix [ll|, |20| 

^ = -\[H,p]-\{Tp + pT)+L{p), (32) 

whereby Zeeman structure of all sublevels of the system 
depicted in Fig. 2 is taken into account. In equation ([32]) . 
the Hamiltonian H describes the system "atom + laser 
field" , the matrix Y describes the spontaneous emission, 
and L{p) describes the cascade effects and has a simple 
explicit form in the representation of polarization mo- 
ments [5J . In order to achieve a fast and efficient solution 
of Eq. (j3"2"|) , we employ the split propagation technique 

mm. 

A. Regular changes of line profiles 

The calculated excitation spectra J{Av) in the case 
of very small laser intensities (Ii as < Icr) are shown in 
Figs. 3(a), 4(a), and 5, whereby the residual Doppler 
broadening has been taken into account by performing 
the convolution (|25|) . An excellent agreement with the 
experimental results is observed. One can also see that 
the relative peak intensities correspond to those expected 
from the theoretical line strengths given in Fig. 6. 

The theoretical spectra in the case when I cr < Ii as < 
I sa t are shown in Figs. 3(b), 4(b), and 5. An interesting 
observation can be made in Fig. 3(b). Intuitively, one 
would expect that optical pumping is manifested more 
strongly and at smaller laser intensities for lines with 
smaller values of branching ratio IT, when only a small 
fraction of population spontaneously returns to the initial 
level. This is, however, not the case. In fact, Eq. (12"3")) 
implies that the critical laser intensity I cr has a minimum 
at II = 0.5. Therefore, nonlinear effects associated with 
optical pumping are more pronounced for HF transitions 
with branching coefficients II close to 0.5. Note, that all 
the excited HF levels considered here have equal lifetimes 

(16.2 ns). Low values of HQ are thus associated with low 

(i) 

values of both the line strengths S{ and the individual 



90 -i — r 
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FIG. 9: (a) The FWHM widths of the lhs (F' = 1) and rhs 
(F' — 2) peaks of Fig. 3 as a function of laser intensity Ii as . 
Solid squares - experiment, F' = 1; open triangles - exper- 
iment, F = 2; solid curves - theory; dashed curves - pure 
power broadening neglecting the broadening due to optical 
pumping. Arrows indicate saturation intensities of both tran- 
sitions, (b) Calculated peak ratio K of the rhs (F' = 2) and 
the lhs (F' = 1) peaks of Fig. 3 as a function of laser intensity 
lias- The calculation was performed with and without taking 
the AmF = ±1 cascades into account. 



Rabi frequencies £1q , such that interaction with laser 
light is inefficient for transitions with small IL\ . Optical 
pumping turns out to be most pronounced for transitions 
with II = 0.5. This can be best seen in Fig. 5: the relative 
intensity of the smaller peak with IT = 1/6 increases 
with respect to the stronger peak with IT = 1/2 as the 
laser intensity is increased in the range I cr < Ii as < I sa t ■ 
This is because the lower level in the case of component 
with IT = 1/2 is faster depleted than it is in the case of 
component with IT = 1/6. 

Another important consequence of optical pumping is 
line broadening, which can be observed when laser in- 
tensity Ii as is close to the critical value I cr [Eq. (23)] 
of the given transition, or when Rabi frequency Q is 
close to the critical Rabi frequency fl cr given by Eq.(21). 
Since I cr has a minimum at II = 0.5, the spectral 
lines with such branching ratio are most strongly af- 
fected by broadening due to optical pumping. Depen- 
dence of the linewidth on laser intensity is illustrated 
in Fig. 9(a) for the 3s 1/2 ,F" = 2 -> 3p 1/2 ,F" = 1 
transition with IT = 5/6 (lhs peak in Fig. 3) and the 
3sy 2 ,F" = 2 -> 3py 2 ,F' = 2 transition with IT = 1/2 
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(a) (b) 




FIG. 10: Zeeman sublevels involved in (a) the 3sx/2,F" = 2 — » 3p 1 /2,F' = 1 transition, and (b) the 3s 1 /2,F" = 2 — > 
3pi/2,F' = 2 transition. Spontaneous emission leads to the population loss to the F" = 1 level of the ground state, and to the 
dark m F n = ±2 levels in (a) and to m F n — in (b). The effective branching coefficient in case (a) is therefore IT e / / (2, 1) < 5/6 
and in case (b) it is II e //(2,2) < 1/2 



(rhs peak in Fig. 3). One can see that remarkable broad- 
ening takes place at laser intensities below the saturation 
intensity I sat [marked in Fig. 9(a) with arrows], while the 
component with 11, = 1/2 (F' = 2) exhibits broadening 
at smaller intensities than the other transition. For com- 
parison, dashed curves in Fig. 9(a) show the intensity de- 
pendence of the power broadened linewidths calculated 

as Av pow — Av nat \J 1 + I 2 /I 3]. It is immediately 
obvious that in the laser intensity range considered here 
the power broadening is much smaller than broadening 
due to optical pumping even at intensities exceeding the 
saturation intensity. 



B. Irregular changes of line profiles and Zeeman 
structure 

Variations of the excitation spectrum of the 
3s 1 / 2 ,F" = 2 — > 3pi/2,F' = 1,2 transitions with 
laser intensity (Fig. 3) are significantly different from 
those observed for the 3si/ 2 ,F" = 1 — > 3pi/2,F' = 1,2 
transition (Fig. 5) in two ways: (i) relative intensity 
of the peak with H = 1/2 first decreases slightly and 
then increases as the laser intensity is increased; (ii) 
Broadening of the peak with II; = 1/2 is actually smaller 
than broadening of the peak with H = 5/6. The key of 
understanding such striking differences is in the different 
Zeeman sublevel structure in both cases. Numerical 
simulations using Eq. (|32p . which include cascade 
transitions with Amp = ±1, yield a ratio 5ft between 
the peak with H = 1/2 and the peak with II,; = 5/6, 
which initially decreases with increasing laser intensity 
and reaches a minimum at Ii as ps 1.3 mW/cm 2 [see 
Fig. 9(b)]. As laser intensity is further increased, the ra- 
tio starts growing, reaches unity at Ii as « 3.2 mW/cm 2 , 
and grows to values slightly larger one. If the cascade 
transitions with Amp = ±1 are ignored, the calculations 



yield a monotonously decreasing ratio 5R [lower curve in 
Fig. 9(a)] without any "abnormalities". 

The effect of Araj? = ±1 transitions becomes obvi- 
ous at closer inspection of Zeeman sublevels involved 
in the 3si /2 ,F" = 2 -► 3p l/2 ,F' = 1 [Fig. 10(a)] and 
is l/2 ,F" = 2 -► 3px/ 2 ,F' = 2 [Fig. 10(b)] transitions. 
Since the laser field is linearly polarized, only the lev- 
els with the same mj? are coupled by it. The presence 
of "dark" levels becomes immediately obvious. In the 
case of the F" = 2 — > F' = 1 transition, the mp = ±2 
sublevels of the lower level are not coupled by the laser 
field and thus act as dark states, which accumulate pop- 
ulation channeled to them vie optical pumping from the 
mp = ±1 sublevels of the upper level. As a result, the 
branching coefficient II; = 5/6 should be replaced by a 
smaller effective branching coefficient H e ff(F" =2,F' = 
1) < 5/6, which accounts for the population loss to dark 
states. In the case of the F" = 2 — > F' = 2 transition the 
dark state is mp = (due to the selection rule Amp ^ 
for F" = F'), therefore II e// (2,2) < 1/2. 

Increase of laser intensity leads to a larger population 
of the dark states, and, consequently, to a monotonous 
decrease of H e ff(F",F'). At very weak laser fields 
II e// (2,l) = 5/6, n e// (2,2) = 1/2, and 7 cr (2,l) > 
Z cr (2,2) [see Eq. ([2"5])]. Hence, the transition F" = 2 -> 
F' = 2 is more strongly affected by optical pumping than 
the other transition. Correspondingly, the ratio 3? de- 
creases with increasing laser intensity. As laser intensity 
is further increased, both II e //(2, 1) and n e //(2, 2) de- 
crease. At some value of Ii as both effective branching 
ratios satisfy the equality II e //(2, 1)(1 — IT e //(2,l)) = 
n e// (2,2)(l - n e// (2,2)). In that case, 7 cr (2, 1) = 
7 cr (2, 2), and both HF transitions are equally strongly 
affected by optical pumping. This corresponds to the ra- 
tio 3? = 1 at Ii as = I eq = 3.2mW/cm 2 in Fig. 9(b). At 
lias > I eq the value of II e //(2,2) becomes larger than 
II e /y(2,l), such that Z cr (2, 1) < Z cr (2, 2), and the ratio 
3? becomes larger than one. At large laser intensities 
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the ratio 3? asymptotically approaches the value of 1.09. 
This is because the populations of the dark mpn levels 
reach their maximum possible values when other mj?» 
levels are fully depleted. In the large intensity limit the 
values II e //(2, 2) and II e //(2, 1) differ by only 9%, there- 
fore both HF transitions exhibit similar broadening due 
to optical pumping (see Fig. 3). In contrast, in the case 
of the 3si/2,F" = 1 — ► 3pi/ 2 ,F' = 1,2 transitions the 
component with IF = 1/2 is apparently more strongly 
broadened than the component with IF = 1/6, which 
could be expected (see Fig. 5). 

VIII. SUMMARY 

We have analyzed the effects of line broadening and 
redistribution of relative peak intensities in the hyper- 
fine excitation spectra of Na atoms due to optical pump- 
ing in the weak excitation limit, when interaction times 
of atoms with the laser field are long compared to the 
characteristic optical pumping time. The study was mo- 
tivated by the lack of availability of detailed theoretical 
models describing such kind of effects in partially open 
level systems at laser intensities below the saturation 
limit. A number of significant results were obtained: (i) 
it is shown that spectral lines can be significantly broad- 
ened at laser intensities well below the saturation inten- 
sity, which is usually regarded as a threshold for onset of 
broadening effects; (ii) it is shown that the presence of 
dark tuf sublevels can vary the effective branching co- 
efficients of the transitions, and this variation depends 
on laser intensity. Changes in the effective branching 
coefficients lead to irregular changes of peak ratios, like 



minimum in the intensity dependence of the peak ra- 
tio, which deviate from those expected from the given 
original branching coefficients; (iii) analytical expressions 
are derived, which allow the calculation of critical val- 
ues for laser intensity and Rabi frequency, above which 
linewidths and peak ratios are notably affected by optical 
pumping; (iv) it is shown that the critical laser intensity 
and critical Rabi frequency depend on the branching co- 
efficient II of the transition, and they have a minimum 
at II = 1/2. 

Accurate theoretical simulations of the density matrix 
equations of motion using the split propagation technique 
yielded a good agreement with the experimental obser- 
vations. In this study we have explored the limiting case 
of long interaction times of atoms with laser field, which 
justified the use of the adiabatic elimination approach. 
It is possible, however, to obtain explicit formulas for 
the excitation spectra in the weak excitation limit also 
without the limitation of adiabaticity in switching Gaus- 
sian laser pulses. In the forthcoming publication we shall 
discuss some unexpected effects related to transit time 
broadening in the other limiting case, when the transit 
time is much smaller than the natural lifetime. 
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